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Abstract 

Two-loop renormalization-group analysis for a two-dimensional electron sys- 
tem with a partially flat Fermi surface has been carried out. We found that, 
irrespective of pairing mechanism, the quasiparticle weight is anisotropically 
suppressed due to logarithmically singular processes in the flat regions of the 
Fermi surface. When the energy scale decreases, the quasiparticle weight is 
the most strongly suppressed around the center of the flat region, which qual- 
itatively agrees with the anisotropic pseudogap behavior suggested through 
the angle-resolved photoemission spectroscopy experiments for underdoped 
Bi2Sr2CaCu208. 
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Non-Fermi-liquid behavior in the normal state of high T c superconductors has provoked 
a great deal of controversy.S Especially, the pseudogap phenomena in the underdoped region 
has been one of the main points at issue. Recent angle-resolved photoemission spectroscopy 
(ARPES) experiments on underdoped B^S^CaC^Og!'! have revealed evidence that on 
cooling the pseudogap determined through the leading edge of the ARPES spectra first 
opens up around the points (±7r, 0) and (0, ±7r) in the Brillouin zone below some crossover 
temperature T* ~ 180K without any broken symmetry. Then, the Fermi surface is broken 
up into the disconnected arcs which shrink with decreasing temperature above the super- 
conducting transition temperature T c ~ 85K. 

Self-consistent treatments of the ci-wave pairing fluctuation above T c showed that the 
Fermi surface is broken up into the patches around (±7r, 0) and (0, ±7r).@ Resonating valence 
bond picture also predicts the patchy Fermi surface due to spinon pairing.! Here we propose 
another mechanism. 

One aspect to study interacting fermion systems is to see how physical quantities flow 
toward their low-energy asymptotic behavior by means of the renormalizat ion-group (RG) 
framework.! Since the ARPES experiments detect the coherent one-particle excitations, 
the key issue here is the quasiparticle weight, z(k), of an electron with a momentum k. 
To see the RG flow of the quasiparticle weight, we have to carry out the RG analysis 
at least up to a "two-loop" order, since the lowest order perturbative correction to the 
quasiparticle weight comes from the logarithmically singular two-loop self energy diagram.! 
In this work, by applying the two-loop RG analysis to a two-dimensional (2D) interacting 
electron system with a partially flat Fermi surface, we present a new scenario, irrespective 
of pairing mechanism, that the quasiparticle weight is anisotropically suppressed due to 
logarithmically singular processes in the flat regions of the Fermi surface. 

We consider a 2D electron system with the model Fermi surface shown in Fig. 1 which 
resembles the real Fermi surface of B^S^CaC^Ogl'! and consists of flat regions (a, a, f3, (3) 
with length 2 A and round-arc regions (A, A, B, B). The centers of the flat regions correspond 
to (±7r, 0) and (0, ±7r) points in Bi 2 Sr 2 CaCu20g. The same model Fermi surface has been 
studied within the one-loop RG analysis (parquet approach) by Zheleznyak, Yakovenko and 
Dzyaloshinskii.i 

Following the standard RG technique,! we consider the linearlized one-particle dispersion 
around the flat regions of the Fermi surface with bandwidth cutoff Eq. We assume Eq <C VpA 
with vp being the correspondign Fermi velocity. As in the case of one-dimensional (ID) 
metals,0 the electron one-particle and two-particle processes within the parallel flat regions 
give rise to the logarithmically singular contribution, ln[-Eb/k>] with uj (< Eq) being an 
infrared cutoff energy scale, to the two-particle scattering strengths and the quasiparticle 
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weight. On the other hand, the processes involving at least one electron in the round- 
arc regions and the processes involving two electrons in two flat regions perpendicular to 
each other, for example a and (3, never give rise to this type of singular contribution. 
Consequently, for an electron located in the round-arc regions, the quasiparticle weight z(k) 
remains marginal. Thus only the processes which are within the parallel flat regions are 
relevant to the RG flow of z(fe).i 

We start with the action for the electron one-particle and two-particle processes in the 
parallel flat regions, a and a, Si = S'km;Z + S'mt;Z, where I is a scaling parameter defined below. 
The kinetic action is given by 

Skin;, = E E r ^ / ^fekXiif)^), (1) 
», Z J-oo 111 JC„., 27T r 



u=a,a a 



where c* CT and the Grassmann variables representing the electron with the spin a 

in the flat region v = a or a and K = (k, e) with e being a Fermion thermal frequency. 
The non-interacting one-particle propagator is given by Q~ l (k x ,ie) = ie — e u (k x ) where 
the linearized one-particle dispersions are given by e a {k x ) = vp{k x — k F ) and Saik x ) = 
vp{—k x — kp) with Vp and kp denoting the corresponding Fermi velocity and the Fermi 
momentum, respectively. We neglect the dependence of the Fermi velocity on k y which is 
irrelevant in the RG sense. We parametrize the infrared cutoff energy as io\ = E e~ l with 
a scaling parameter, I. Two-dimensional electron momenta, k, are restricted to the set 
C v . x = {k || e v (k x ) \< -A < k y < A}. 

The action for the two-particle processes are written as 

Sint-,1 = n~vp E II / |r I 77^ ( 2 ) 

W i=hiJ-°° b ( 2?r ) 2 ~ 2^ -I- 00 2^ h** (2tt) 2 Jc a[l (2tt) 2 ic a;i (2tt) 2 

9{uii/2— |ea(fca!4)|) 9 (A— \k yi \) with fc 4 = ki + k 2 — k 3 . The two-particle scattering processes 

considered here are depicted in Fig. 2. Each leg of the scattering vertex is labeled by the 

regional index v — a, a, the momentum in the k y direction and the spin. The dimensionless 

scattering strengths are given by 

9\ {^VV ^2/2' ^2/3) — ^0-40-2 ^0-30-1 9l (kyii ^V2i ^2/3) ~~ ^0-401 S a2C r 3 Q\ {ky\ i k y2 i ^3/3) (3) 

where g\ l \k yv k y2 , k y3 ) and gf\k yv k V2 , k y3 ) denote the strengths of the backward and for- 
ward scattering between the parallel flat regions, respectively. Within the present model, the 
electron umklapp process occurs in the special case where the distance between the parallel 
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flat regions is equal to n. In this work, we do not take account of this case and neglect the 
electron umklapp process. 

We split up the set of fc-points, C v -i, into two subsets as C v -i = C< l+dl © dC>. l+dl , 
where C< l+dl = {k \\ e u (k x ) \< ui +dl /2, -A < k y < A} (low-energy shell) and dC> l+dl = 
{k | uoi +d i/2 <\ e u (k x ) \< uji/2, —A < k y < A} (high-energy shell). Accordingly, the action 
is decomposed as Si = Sf +dl + Sf +dl . Integration over the modes in the high-energy shell 
gives the corresponding counterpart of the partition function as 



L 



Vc*Vcexp 



oo 

$l+dl + X rr((['5'int;/+^] 



n=l n! 



(4) 



where Vc*Vc symbolizes the measure over the Fermion Grassmann variables and 
C< +dl means that Fermion momenta are restricted to the low-energy shell C< l+dl or 
C^. l+dl . The average over the modes in the high-energy shell is defined as (((• • •))) = 

Z- 1 Vc*Vcexp[S^ in . l+dl ] (•••), with Z> = / Vc*Vcexp[S^ in . l+dl ] and the sub- 

JdC i+di JdC i+di 
script 'c' represents the connected diagrams. We perform a perturbative expansion of 

X^Li ^r(([ , S'int /+(iJ ra ))c by picking up the Feynmann diagrams whose contribution is in pro- 
portion to dl and then replacing S^ +dl + X^Li ^ {{[S^i+diY 1 )) c with the renormalized action 
generally written in a from 

/°° dp c dk 
, - -°o Z7T JC u . l+il [Z7T) 

+™f e n f §L ^[9i im (k yv k y2 , ky3 ) (5) 

+W p^(k yv k y2 , k y3 )dl + 0(dl 2 )} c^K^iK^i^c^i^). 



The contribution Oi(k y )Q v 1 (k x , ie)dl in the renormalized kinetic action comes from the two- 
loop self energy diagram depicted in Fig. 3, where 9i(k y ) is given by 

9^) = I [2gVngU + 2g?>Ug\*> - g^Ug^ - g^Ug™] , (6) 
where the operation represented by □ is defined as 

9? ] Vg\ j) = ^ r M J M %?{k' y , k y , k y - qy )g?{k' y + q y , k y - q y , k y ). (7) 

J-A+k y Z7T >/-A-Min(0,q y ) 271 

The integration ranges are determined through the condition that all the incoming and 
outgoing momenta at each scattering vertex in Fig. 3 must be within the parallel flat 
regions. In the renormalized action for the two-particle scattering processes, the con- 
tribution Wi '^ \k y ^-, k y2 , ky^dl, where Wi 1 2 3 4 (k y ^, k y2 , k y3 ) = 5 cri(72 5i 7a i 7l Wi ^ (k y ^, k y2 , k y ^) 

(2) 

— 5cr 2 ct 3 u>i (k yv k V2 , k y3 ), comes from the two-loop vertex correction diagrams depicted 
in Fig. 4. As space is limited, we do not write down full expression for w\ l \k yi , k yv k y3 ). 
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To restore the original cutoff, we rescale the momenta and frequencies as k = (e dl k x , k y ) = 
(k x , k y ) and e = e dl e, respectively. Then, to leave the kinetic action scale- invariant, we 
simultaneously perform the scale transformation of the Fermion Grassmann variable as 
c ua (K) = [1 + \{0i(k y ) — 3}dl)c va (K), where K = (k x ,k y ,ie). Thus the two-loop RG 
equations for the scattering strengths are given by 



d"9l (ky-^i ky2i ^1/3) (i) , , , l, \ ^ W/i, U U \ \^ 0(U \ (Q\ 

^ — w i \ K yv K y2i K y?,> ^ y^ yi' yv 2^ u A K yj)i \°) 



where i = 1, 2 and k yi = k y \ + k y2 — k y3 . The one- loop counterparts of the RG eqs. (|S]) are 
given by Zheleznyak, Yakovenko and Dzyaloshinskii,i where w\ l \k yi , k V2 , k y3 ) includes only 
one-loop vertex correction diagrams, Fig. 4(a) and 4(b), and the second term of the r.h.s 
of the RG eq. (§) does not appear. As an initial condition at / = 0, we put the scattering 
strengths equal to the momentum-independent Hubbard repulsion, U, i.e. 

9o\k yv k V2 , k y3 ) = U/nv F . (9) 

The RG equation for the quasiparticle weight is obtained through the renormalized ki- 
netic action in eq. (El) and is written as 

d\nzi(k y ) 



dl 



Mky). (10) 



As an initial condition at / = 0, we put the quasiparticle weight equal to the non-interacting 
value, i.e. 

zo(ky) = 1. (11) 

It is here instructive to see the k y dependence of zi(k y ) by assuming that gf\ky V k y2 ,k y3 ) 
remain at the initial value, U/ttvf. This assumption applies to the early-stage (i.e, high- 
energy) behavior of the RG flow. Then, eqs. (H), (0) and ( |T0D give the result, zi{k y ) = 
[u;;/-Eo] c/2( ' 3A2_fc ^^ 167r2 with U = U/ttvf, which indicates that the quasiparticle weight at the 
center of the flat region (k y = 0) is the most strongly suppressed due solely to the kinematical 
restriction to the logarithmically divergent processes. In reality, however, we cannot obtain 
the 2- loop RG flow of zi(k y ) in a consistent manner, until we take account of the anisotropic 
RG flow of the scattering strengths. 

We solve numerically the RG equations, (|8|) and fllCf), under the initial conditions, (0) 
and (|TT]) . Absolute value of the Hubbard repulsion, U, is not essential for the RG flow of 
the scattering strengths and the quasiparticle weight. From now on, we set the Hubbard 
repulsion at U = 4:.5tivf. Linear differential equations, (|) and (|TUD, are treated as recursion 
equations with an infinitesimal difference step dl = 0.08. The numerical procedure at the 
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n-th step with the corresponding scaling parameter l n = ndl consists of three consecutive 
steps. First, we insert g ( f l {k yv k y2 , k V3 ) into the expressions of Oi n {k y ) and w^(k yi , k y2 , k V3 ). 
Second, we perform the integration over the momenta in the k y direction by dividing the 
interval, —A < k y < A, into 32 discrete points to obtain 6i n (k y ) and (k yi , k V2 , k V3 ). 
Finally, we obtain gij +1 (k yv k V2 , k V3 ) by using eq. (|). 

In Fig. 5(a), we show the RG flow of the quasiparticle weight, Z[(k y ), for < / < 3.2. 
When the energy scale decreases, the quasiparticle weight is the most strongly suppressed 
around the center of the flat region, k y = 0, forming a dip there. In Fig. 5(b), we show the 
same RG flow for 3.2 < / < 4.08. The k y dependence of the quasiparticle weight undergoes 
a crossover behavior toward its low-energy asymptotics. There is some crossover scaling 
parameter Z cross ~ 3 around which the dip of zi(k y ) around k y = becomes flatter. Then, 
in the low-energy regime, / > Z C ross, the region around which Zi{k y ) is the most strongly 
suppressed moves from the center [k y = 0) toward the edges {k y = ±A) and finally zi{k y ) 
approaches zero everywhere for —A < k y < A in the low-energy limit. The k y - dependence of 
the quasiparticle weight in the high-energy regime corresponds to the anisotropic pseudogap 
behavior in underdoped Bi 2 Sr 2 CaCu208 suggested by the ARPES experiments. In the 
present analysis, the crossover value of the scaling parameter, / croS s, corresponds to the tem- 
perature scale, T cross ~ 0.05E . Thus we may expect that the temperature range coverved by 
the ARPES experiments on the normal state of Bi 2 Sr 2 CaCu208 lies within the high-energy 
regime introduced here. We here again stress that, in the present scheme, the anisotropic RG 
flow of the quasipartcile weight occurs due solely to the logarithmically singular processes 
coming from the flat regions of the Fermi surface. 

To see the reason why the crossover behavior of zi(k y ) occurs, in Fig. 6, we show the RG 
flow of the backward scattering strength, g\ l \k yi , k y2 , k V3 ) at points, (k vv k y2 , k V3 ) = (0, 0, 0), 
(A, — A, A), (A, — A, — A), (A, A, A). As an example, let us see the RG flow of ^ (1) (0, 0, 0). The 
monotonic decrease of g^\o, 0, 0) at the early stage of the renormalization reminds us of the 
RG flow of the backward scattering strength in the purely one- dimensional metal, where the 
backward scattering strength monotonously decreases and vanishes toward the low-energy 
asymptotics corresponding to the Tomonaga-Luttinger fixed point. I In the present case in 
contrast to this, the RG flow of gj'^O, 0, 0) has a very shallow minimum around I ~ I* and 
begins to increase for I > / croS s- As is seen in the inset in Fig. 6, the flow becomes nearly 
constant for I > l c ~ 4.1 in the low-energy limit due to the two- loop vertex corrections, while, 
in the case of the one-loop RG analysis, the corresponding flow diverges around I ~ Z c .i 

The crossover behavior of g\ l) (0,0,0) originates from breakdown of the cancella- 
tion between the particle-particle (Cooper) loop [Fig. 4(a)] and the particle- hole (Peierls) 
loop[Fig. 4(b)]. As the renormalization process goes on, the breakdown becomes more and 
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more conspicuous and the Peierls loop tends to dominate the Cooper loop. Consequently, 
the spin-density-wave fluctuation becomes the most dominant fluctuation in the low-energy 
limit, provided that the regions a and a are perfectly flat.! The crossover behavior of the 
RG flow of zi(ky) corresponds to the crossover behavior of the scattering strengths from the 
high-energy regime to the low-energy regime. 

We here consider the effects of small but finite curvature of the Fermi surface in the 
regions a and a which always exists in reality. In this case, at the energy scale smaller 
than the energy resolution, AE, which detects the Fermi surface curvature, the two-loop 
self-energy diagram never gives rise to the contribution of the type Oi(k y )Q~ 1 (k x , ie)dl in the 
renormalized kinetic action in eq.fl^). Consequently, the RG flow of zi(k y ) stops around the 
scaling parameter, l cuvv , defined as AE ~ E e~ lcuiv . In Fig. 7, we take Z curv = 3 (correspond- 
ing to AE/ Eq ~ 0.05) as an example and show the RG flows of Zi(k y ) at some representative 
points on the Fermi surface. As the energy scale (or temperature scale) decreases, the quasi- 
particle weight at the point a on the round-arc region remains marginal, indicating that 
the round-arc regions of the Fermi surface are robust. In the flat regions, the quasiparticle 
weight at the center, c, is more strongly suppressed than that at the edge, b, but remains 
finite in the low-energy limit due to the curvature of the Fermi surface. 

In our analysis, the RG flow of Zi(k y ) becomes discontinuous at the edge points of the 
flat region (the point b in the inset of the Fig. 7), which seems unnatural. We give here a 
qualitative comment on the origin of this discontinuity. In the present analysis, we treated 
the flat regions as decoupled to the round-arc regions in the RG sense. However, strictly 
speaking, at finite energy scale, the flat regions and the round-arc regions can couple to 
each other through the scattering processes involving the very tiny region, the shaded area 
in Fig 8, which can be regarded as an extension of the flat region. Thus, if we take account 
of these processes, the discontinuity of zi(k y ) at the edge points will be removed. We do 
not inquire further into this point, since this superficial discontinuity is irrelevant to our 
qualitative understanding of the RG flow of zi(k y ) in the present model. 

In summary, we have carried out the two-loop renormalization-group analysis for the 
two-dimensional electron system with the partially flat Fermi surface and found that the 
logarithmically singular processes coming from the flat regions cause the renormalization- 
group flow of the quasiparticle weight which depends on the location of the fc-point in the 
flat regions of the Fermi surface. When the energy scale decreases, the quasiparticle weight 
is the most strongly suppressed around the center of the flat region, forming a dip there 
which qualitatively agrees with the anisotropic pseudogap behavior suggested through the 
ARPES experiments. 
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FIGURES 



FIG. 1. The model Fermi surface considered here. 

FIG. 2. The two-particle scattering vertex considered here. The solid and broken lines represent 
the propagators for the electrons in the flat regions, a and a, respectively. The box represents the 
scattering vertex, g^ ia2a3<7 * (k yi , k y2 , k y3 ). 

FIG. 3. The two-loop self-energy diagram. 

FIG. 4. The diagrams which contribute, (a) to the one-loop vertex correction in the parti- 
cle-particle (Cooper) channel, (b) to the one-loop vertex correction in the particle-hole (Peierls) 
channel, and (c) to the two-loop vertex corrections. 

FIG. 5. k y dependence of the renormalization-group flow of the quasiparticle weight, zi(k y ), 
(a) for < / < 3.2, and (b) for 3.2 < / < 4.08. 

FIG. 6. The renormalization-group flow of the backward scattering strength, g^ 1 \k yi ,k y2 , k y3 ) 
at (k yv k y2 , k y3 ) = (0,0,0), (A, —A, A), (A, —A, —A), (A, A, A). In the inset, the same graph 
covering a wider range of the vertical axis is shown. 

FIG. 7. The renormalization-group flow of the quasiparticle weight, zi(k y ), at some represen- 
tative points on the Fermi surface. 

FIG. 8. The shaded area also contributes to the renormalization-group flow of the quasiparticle 
weight, which is not taken into account in this study. 
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